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ABSTRACT 

We investigate the effect of modified gravity on cluster abundance and the Sunyaev- 
Zel'dovich angular power spectrum. Our modified gravity is based on a phenomeno- 
logical extension of the Dvali-Gabadadze-Porrati model which includes two free pa- 
rameters characterizing deviation from ACDM cosmology. Assuming that Birkhoff's 
theorem gives a reasonable approximation, we study the spherical collapse model of 
structure formation and show that while the growth function changes to some extent, 
modified gravity gives rise to no significant change in the linear density contrast at 
collapse time. The growth function is enhanced in the so called normal branch, while in 
the "self-accelerating" branch it is suppressed. The Sunyaev-Zel'dovich angular power 
spectrum is computed in the normal branch, which allows us to put observational 
constraints on the parameters of the modified gravity model using small scale CMB 
observation data. 
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1 INTRODUCTION 

General Relativity is surely the most successful theory of 
gravity that passes accurate tests in the solar system and 
laboratories. However, current cosmological observations in- 
dicate the presence of dark matter and dark energy; a large 
fraction of the Universe is made of unknown components. 
The mystery of the dark components is based on general 
relativity, and hence it tells us that what we do not know 
may be the long distance behaviour of gravity rather than 
the energy-momentum components in the Universe. In this 
sense, cosmological observations open up a new window to 
study the properties of gravity on large scales. 

There are various alternative theories of gravity lead- 
ing to interesting cosmological consequences. For exam- 
ple, t he Dvali-Gabadadze-Porrati (DGP) model (|Dvali et al.l 
2000) is one of the extra dimensional scenarios that can ac- 
count for cosmic acceleration without introducing dark en- 
ergy. The f(R) theories, which modify the four-dimensional 
Einstein-Hilbert action explicitly, also reali se the ac- 
celera ted expan sion of the U niverse (|Sotiriou fc Faraonil 
|200S|). MOND (iMilgroml If 983T ) and its relativistic exten- 
sion |Bekensteinll2004 ) explain galactic rotation curves with- 
out need for dark matter. A more phenomenological way 
of changing gravity is to assume Yukawa-like modification 
to a gravitational potential. Such modification yields effec- 



tively a scale-dependent Newtonian constant, and its ef- 
fect on the evolution of large s c ale structure has b e en in- 
vestigated ilSealfon et all 120051; Ishirata et af] 120051 . 120071 ; 
IStabenau fc Jain 2006; Marti no et"afll2008h . 

One of the powerful ways for distinguishing modified 
gravity from the ACDM model is to study the growth func- 
tion because the modified growth function would leave its 
footprints on the cosmological large-scale structure. In this 
context, many authors have studied the integrated Sachs- 
Wolfe effect and weak lensing in modified gravity, and ob- 
tained obser vational constrai nts on modified gravity models. 
For instance. ISchmidtl l|2008l ) investigated the effect of f{R) 
gravity, the DGP model, and tensor-vector-scalar theory on 
weak lensing, and showed that for detecting signatures of 
modified gravity the weak lensing observation is a better 
probe than the integrated Sachs- Wolfe effect measured via 
the galaxy-CMB cross-correlation. I Thomas et al.l l|2008l ) put 
constraints on the DGP model by using weak lensing data 
(CFHTLS-wide) com bined with baryon ac oustic oscillations 
and supernovae data. ISchmidt et all (2009) have studied the 
statistical properties of dark halos in f(R) gravity by em- 
ploying numerical simulations. 

In this paper, we consider a generalisation of the 
DGP model, adding a term ±_ff 2c */rc' 1_Q ' in the Fried- 
mann equation, where a and r c are the model parame- 
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ters dDvali fc Turned \200A IKovamal 120061 ; lAfshordi et~all 
|200S| ; iKhourv fc Wvmanll2009l ). The generalised model re- 
duces to the self-accelerating and normal branches of the 
DGP model (with a cosmological constant) by taking a — > 
1/2, and reproduces the ACDM model in the a — > 
limit. Af shordi et alJ (|2008h studied cosmological perturba- 
tions in this modified gravity theory under several assump- 
tions, and discussed various ob s ervati onal consequences of 
the model. IKhourv fc Wvmanl l|2009l ) performed N-body 
simulations in the same model recently. We shall work on 
a different approach in the present paper: assuming that 
Birkhoff ' s theorem gives a reasonable approximation (as im- 
plied by IKovamal (|2006h ). we study the spherical collapse 
model of structure formation with the modified Friedmann 
equation. Using the spherical collapse scenario, we can com- 
pute the Sunyaev-Zel'dovich (SZ) angul ar power spectrum 
through the Press-Schechter formalism (|Press fc Schechterl 
Il974l ). Since the SZ angular power spectrum is sensitive to 
the distribution of dark halos, it is a good probe of modified 
distribution of dark halos in the generalised DGP model. 

Before closing the introduction we should note that in 
any case modified gravity must satisfy solar system and lab- 
oratory tests, and some f(R) models and the original DGP 
model indeed have mechanisms to reproduce ordinary grav- 
ity around the Sun and on the Earth. This point is not clear 
in the present model and is beyond the scope of the paper. 
We only remark that our Friedmann equation reduces to the 
ordinary one at high densities. 

This paper is organized as follows. In the next section, 
we describe our modified gravity model in terms of the Fried- 
mann equation and explain its possible origin. In Sec. [3] we 
review the spherical collapse model of structure formation 
in modified gravity, and compute the linear growth func- 
tion and the linear density contrast at collapse time. Then, 
in Sec. |4j we study the effect of modified gravity on the 
SZ angular power spectrum. Observational constraints on 
the model parameters are discussed. Finally, we conclude in 
Sec. 



2 THE MODEL 

The DGP braneworld was originally proposed as a model 
for recovering 4D gravity o n the brane even in an infinitely 
large 5D Minkowski bulk ijDvali et al.l l2000h . In the DGP 
braneworld, gravity is modified at long distances, while 
standard 4D gravity is indeed reproduced at short dis- 
tances (th£ough_aj^ompljrat^ 

Tanakall2004 iKovama fc Silvall2007t iLue fc Starkmani 2003; 



Lue et all 2004 ). so that the model passes solar system and 
laboratory tests. Probably the most intriguing consequence 
of the DGP br aneworld comes out in the modified Fried- 
mann equation (iDeffavet et alj|2002t ): 

H 2 = ^p±?L. (1) 

3 r c 

Here, r c is the crossover scale above which gravity looks 5D. 
With the upper (plus) sign and r c ~ Hg 1 (the present Hub- 
ble horizon scale), we have the "self-accelerating" solution 
which could be the origin of the current cosmic acceleration. 
In this paper, however, we do not assume that modified grav- 
ity is directly responsible for the accelerated expansion. 



The model we consider is a phenomenological extension 
of the DGP braneworld described by the modified Fried- 
mann equation 



H 2 



a 8ttG A 



a) '■ 



w here < a < 1. This is similar to the model 

of iDvali fc Turner! (|2003h . but we allow for a different sign 
of the last term and include t he cosmological constant A 
explicitly (jAfshordi et all l2008h . With the upper (respec- 
tively lower) choice of sign we use the terminology the self- 
accelerating (respectively normal) branch, though it is not 
se//-accelerating in the upper sign case. Equation (J2J with 
q = 1/2 corresponds to the DGP cosmology (with a cosmo- 
logical constant or the tension on the brane), while a = 1 
can be absorbed into a redefinition of the gravitational con- 
stant G. Expansion history of ACDM cosmology is recovered 
in the limit a — ► 0. 

The modified Friedmann equation J2J can be recast in 



JT = — +^(r c H ) 



-2(1-0,) 



where 



A := 1 - fi m o =F (rcHo) 



-2(l-c) 



(3) 



(4) 



fimo := 8-kGpo/(SHq), and the present scale factor is chosen 
to be ao = 1. If A(oc A) = and the accelerated expansion 
were supported by modified gravity in the self-accelerating 
branch, successful cosmology would require r c ~ Hq . In the 
present case, however, A is responsible for the accelerated 
expansion as in conventional cosmology, and therefore we are 
in principle allowed to take r c < Hq 1 (e.g., r c ~ O.IJYq -1 ) [3 
In particular, for sufficiently small a, expansion history is 
very close to that in stan dard ACDM cosmolo gy even if we 
take relatively small r c ( Afs hordi et al.l [2008T ). Indeed, for 
a = 0.1, r c Ho = 0.4, and Sl m o = 0.26, the dimensionless 
physical distance 



E{a) := H [ - 



da 



H( a y 



(•>) 



differs from the corresponding ACDM result by less than 
one percent (Fig. [TJ). 

Unfortunately, we do not have concrete higher dimen- 
sional models that account for Eq. ((2j). One possibility is that 
such modification could be derived fro m a higher codimen- 
sion DGP model (jAfshordi et al.|[200 8). as explained below. 
Let us consider a graviton propagator which is proportional 
to 



(6) 



2(1-7) ' 



k 2 + k 2 -< /r, 

This follows from a phenomen ological model of modified 
gravity proposed in lDvalil (|2006l ) and is a power-law general- 
isation of the graviton propagator in the DGP braneworld. 
Since we are interested in long distance modification of grav- 
ity, we assume tha t 7 < 1. The unitarity constraint requires 
7 > (Dvali 20061). It can be seen that 7 = 1/2 reproduces 



1 The normal branch does not admit acceleration without A from 
the beginning. 
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Figure 1. Dimcnsionlcss physical distance. Parameters are (from 
top to bottom): (1) a = 0.1, r c Ho = 0.4; (2) a = 0.01, r c Ho = 
0.2; (3) ACDM; (4) a = 0.01, r c H = 0.2; (5) a = 0.1, r c H = 
0.4. The cases (1) and (2) are for the normal branch and the cases 
(4) and (5) are for the self-accelerating branch. In all the cases 
n m0 = 0.26. 



the original DGP model and 7 = 1 can be absorbed into a re- 
definition of G. From this observation, we may simply iden- 
tify 7 = a. The propagator with 7 1 has s ome connection 
to the so called cascading DGP braneworld (|de Rham et alj 
l2008al lbl)FI and so the Friedmann equation with a < 1 might 
be realised in such higher codimensi on models (see also 
iKaloper fc Kile y 2007; Kobayashi 20081 ). However, we would 
like to stress that detailed analysis of higher codimension 
DGP models has yet to be undertaken and no braneworld 
models have been known so far that lead to Eq. ||SJ. There- 
fore, we shall view Eq. © as a phenomenological starting 
point of our modified gravity. 



3 STRUCTURE FORMATION IN MODIFIED 
GRAVITY 

In the original DGP case, we have the covariant gravita- 
tional field equations that not only lead to the Friedmann 
equation {TJ but also govern the behaviour of cosmological 
perturbations. As we do not know a complete set of field 
equations that underlies our phenomenologically extended 
model, we must assume something about the dynamics o f 
cosmic inhomogeneities in the present case. iKovamal (|2006l ) 
clarified this issue by constructing simple covariant grav- 
itational equations which give essentially the same Fried- 
mann equation as Eq. ([2]). The gravitational equations con- 
tain an additional term called E^ v in order to satisfy the 
Bianchi identity, which hinders to get closed form equa- 
tions. In the original DGP braneworld, the evolution of the 
E^v term follows from the full 5D Einstein equations. In 
the absence of underlying theories for general a, one has 
to assume the struct ure of E u u , whi ch then determines the 
growth of structure. IKovamal (|2006l ) considered two possi- 
bilities: (i) weak gravity is described by the scalar-tensor 
theory, as in the original DGP model; (ii) the modified 



2 The codimension-two case in fact gives the propagator ~ In k 
for small fc. 



Figure 2. The growth function divided by the scale factor for 
different parameters. Plots are for the normal branch and f2 m o = 
0.26. 



gravity m odel respects B ir khoff 's theorem (at least approx- 
imately). lAfshordi et al.l (|2008l) invoke the par ameterized 
post-Friedmann framework (|Hu fc Sawi cki 2007) and hence 
effectively take the first approach. In this paper, we employ 
the second approach and study nonlinear structure forma- 
tion in modified gravity. The assumed Birkhoff's theorem 
allows us to use the modified Friedmann equation to track 
the nonlinear dynamics in a simple way, rather than intro- 
duce the extra scalar degree of freedom explicitly. The modi- 
fied Friedmann equation reduces to the usual one in the high 
density regime, and in this sense we implem ent the nonlinea r 
recovery of GR. In the cases investigated bv lKovamal l|2006l ). 
the difference between the above two approach is small con- 
cerning the linear growth of perturbations[f| We come back 
to this issue in Appendix. 

In this section, we compute the linear growth func- 
tion and the linear density contrast for spherical collapse, 
whic h are the key quantities for the Press-Schechter formal- 
ism (|Press fc Schechterlll974l ) to predict the number density 
of clusters. 



3.1 The growth equation 

Following Sch aefer fc Kovamal l|2008l ). we study a spherical 
overdensity with matter density p c = p c (t) and radius R = 
R(t) in a background governed by the modified Friedmann 
equation @, which can be recast in 

8ttGp 



H J = H$g(t), £:= 



3H$ 



Differentiation with respect to t leads to 
3, 



m 



9(0 - Ug'(0 



(7) 



(8) 



where a := da/dt and g := dg/d£. 



3 lAfshordi et al ] [|2008h focuses on the a — > limit, but in the 
PPF framework cosmological perturbations are still sensitive to 
r c . On the other hand, the Birkhoff's theorem-based approach 
relics essentially on the Friedmann equation, and hence the effect 
of modified gravity vanishes in the a — > limit. Thus, the two 
approaches give different predictions at least in this limit. 
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Figure 3. The growth function divided by the scale factor for 
different parameters. Plots are for the self-accelerating branch 
and S7 m o = 0.26. 



Our central assumption is that the modified gravity the- 
ory respects Birkhoff's theorem. The dynamics of R(t) is 
then described by 



R 
R 

with 



Hi 



8nGp c 



pc oc r 



We define the overdensity as 

S ._ Pc - P 



(9) 



(10) 



(11) 



Upon linearisation the evolution equation for 8 is given by 

6 + 2H5 = 4nG [g'(0 + 3£g"($] pS. (12) 

The growth function D+(a) is defined by <5(a,x) = 
D+(a)S(ao, x). It follows from Eq. (JT2J) that 

da^ a da 



D 4 



where 



fi = 3 + 



dlnH 

din a ' 



dlntf /din// 

/2 = O— h 



d In a 



d In a 



+ 



a 2 d 2 // 
//da 2 "' 



(13) 

(14) 
(15) 



The boundary condition is given by D+(0) = and 
D+(ao) = 1. Note that in standard cold dark matter domi- 
nated cosmology we find D + (a) = a. 

Figures [2] and [3] show typical behaviour of the growth 
function. For example, the difference between background 
expansion history E(a) in modified gravity with (a, r c Ho) = 
(0.01, 0.2) and that in the ACDM model is < 3 %, for which 
the growth functions differ by almost 10 % (both in the 
self-accelerating and normal branches). In the case of the 
self-accelerating branch, the growth of perturbations is sup- 
pressed comp ared to the ACDM model, in agreement with 
the result of ISchaefer fc Kovamal (120081 ) H This is because 
the effective dark energy term becomes dominant earlier 

4 Note that the normalization of the growth function here is such 
that D-|-(ao) = 1. Hence, smaller (respectively larger) D+(a) im- 



Figure 4. Overdensity at collapse time for different parameters. 
Plots are for the normal branch and f2 m o = 0.26. 
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Figure 5. Overdensity at collapse time for different parameters. 
Plots are for the self-accelerating branch and Q m o = 0.26. 



in modified gravity than in the ACDM model due to the 
+H 2a term. In the normal branch, the effect of modified 
gravity works oppositely and the growth of perturbations 
is enhanced compared to the ACDM model. In the both 
cases, deviation from the ACDM result is larger for larger 
a and smaller r c , as expected. Note here that the PPF ap- 
proach predicts qualitatively the same result: structure is 
more evolved in the norm al branch than in ACDM cosmol- 
ogy (|Afshordi et al.|[2008t ). Note also that qualitatively the 
same result, i.e., the suppressed (enhanced) growth func- 
tion, is also found in the self-accelerating (normal) branch 
of the original DGP model (a = 1/2) by solving th e five- 
dime nsional Einstein equations (| Cardoso et al. 2008; ISond 
l200ct) . 



3.2 Spherical collapse 

In order to study spherical collapse, it is convenient to use 
the quantities normalised by their values at turn-around 



time dWang fc Steinhardtlll998l ; iMota fc van de Bruckll2004l ; 
iBartelmann et all [2006). First, we define the normalised 



plies a smaller (respectively larger) amplitude of the perturbation 
at a < 1 evolving into some fixed amplitude at a = 1, i.e., the 
growth is enhanced (suppressed). 
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scale factor and radius of the overdensity as follows: 

s:=a/a ta , y := R/Rta.. (16) 

We also define the dimensionless time r := //tai, where 
Hta, := H(ata) is the Hubble rate at turn around. Now the 
modified Friedmann equation can be written as 

2 

= X + 1 - w 



HI 



±{Hujr c ) 



-2(l- a ) 



(17) 



where lo := Q.^Hq j a 3 a H 2 ^, \ '-~ u/x 3 , and a dot here and 
hereafter denotes derivative with respect to r. This equation 
can be rewritten as (x/x) 2 = h{x), or, equivalently, 



x — \ I x 2 h \ —? 



(18) 



Similarly to the previous calculation, one obtains the an- 
other Friedmann equation that describes the evolution of 
the overdensity patch: 



y = y 



3 ( 



2 </3 



V 



(19) 



where b! := dh/d\ and ( := (p c /p) \x-i- The boundary con- 
dition y\ x =o = 0, y\ x =i = 0, and y\ x =i = 1 uniquely deter- 
mines Equation (I19p reduces to a first order differential 
equation by noticing that 

d . 2 _d_ 

dr 



d7 y 



2, / 



From this we obtain 

tf = y 2 h(^f)-K&), (20) 

where we fixed the integration constant by using the bound- 
ary condition y — at turn-around time (y = 1). 

Since the background dynamics at early times is the 
same as that of the standard matter dominant universe, we 
may approximate h(w/x 3 ) ~ lo/x 3 for i < 1. Thus, at early 
times we simply have 

Ix 3 ' 2 



Similarly, Eq 
leading to 

. 2 v m 



20)1 reduces to 

KC u ) y 



3 a/Cw 



1 + 



(lj 2 



3 h(Cu ) 
10 



dj/ for y <§C 1, 



(21) 



(22) 



(23) 



The time evolution of the nonlinear overdensity, A := 
£x 3 /y 3 , can be computed at early times by using 



1+ 3^) 
5 Qlj 



(24) 



The linear density contrast at collapse time, 8 C , can be used 
to relate the nonlinear overdensity A with the density that 
would result from linear evolution and the same initial con- 
dition. Using the growth function, we have 

D+(x c 



lim , 

x^o D+(x) 



[A(x) - 1] . 



(25) 



We numerically computed 5 c (z) for various model pa- 
rameters, and the results are illustrated in Figs. [4] and [5] 
We find that modified gravity does not change 8 C much: for 
example, S c in modified gravity with (a, r c Ho) = (0.01, 0.2) 
differs from the ACDM prediction by 1 - 2 % (both in the 
self-accelerating and normal branches) 



4 THE SUNYAEV-ZEL'DOVICH POWER 
SPECTRUM 

In this section, we investigate the effect of modified gravity 
on the SZ angular power spectrum. Modification of gravity 
changes cluster number count through the modified growth 
function and critical density contrast. We concentrate on the 
case of the normal branch. The SZ angular power spectrum 
could be amplified in this branch because the growth func- 
tion is enhanced relative to the ACDM model. The change 
in the critical density contrast is quite small, and hence it 
will give rise to only a negligible effect. 

Modified gravity might in general affect the halo 
profiles assumed in the following discussion. This is- 
sue can be addressed by evaluating the Vainshtein ra- 
dius, r„, below which Einstein gravity is recovered. The 
Vainshtein radius in this particular model is given by 
r, = [ r *( 1 - a ) rs ] 1 /[i+4(i-a)] ~ [ r 4 r g]V5 w here r s is the 
Schwarzschild radius of the source (|Dvalill2006h . More ex- 
plicitly, one has r, ~ W 2 [(r c H ) 4 (M/M e )} 1/5 kpc, and for 
r c H ~ 0.01 and M ~ 10 14 M Q , r, ~ IMpc. This estimate 
allows us to employ the halo profiles used in the context of 
standard gravity. 

The computation of the SZ angular power spec- 
trum is based on the ha l o formalism |Cole Sz Kaiser] 
19881: iMakino fc Sutol Il993l : iKomatsu fc Kitavamal Il999l : 
Komatsu fc Sel iak 2002). Since the one-halo Poisson term 
dominates the halo-halo correlation term on the scales we 
are interested in, we neglect the halo-halo correlation term. 
The SZ angular power spectrum is then given by 



r sz _ 2 
1st — 9v 



dz 



dz 



^M^fl M M, Z )\ 2 ,(26) 



where V(z) is the comoving volume at z per steradian, 
dn(M, z)/dM is the number density of clusters, yi(M,z) 
is the 2D Fourier transform of the projected Compton y- 
parameter, and (?„ is the spectral function of the SZ effect, 
which is given by 



9v = 



I) 2 



a;coth ( —x 



(27) 



where x = hv /fcsT 7 with the Boltzmann constant fes and 
the CMB temperature T 7 . 

As we are considering the spherical collapse scenario 



5 It was reported that in the original DGP model with A = 0, 
the overdensity for spherical collapse decreases by almos t 20 % 
compared to conventional cosmology (Schacfer & Koyama 2008). 
However, the first version of the preprint (i.e., arXiv:071I.3I29 
vl) includes an error in the numerical calculations. We confirmed 
with the authors that <5 C differs from the s tandard one only by 
a sma ll amount also in the case studied by ISchaefer fc Kovamal 
(120081) . 
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for the formation of halos, we utilise the Press-Schechter 
theory (|Press fc Schechtedll974h , 



dn(M, z) 
dM 



ttM 



do 



t{M,z) dM 





" &l 1 


exp 


2a(M,z)_ 



(28) 



where 5 C is the critical over density which is obtained in the 
previous section, a(M, z) is the variance of the matter den- 
sity field on the mass scale M. The variance a is computed 
from the power spectrum of linear matter density fluctua- 
tions with the top hat filter, 



(29) 



a(M,z) = / &kk z P m {k,z)W{kR), 



where W(kR) is the top hat window function and R is the 
scale which corresponds to M(— AnpR 3 /3). We calculate 
the linear power spectrum P ul (k,z) in the modified gravity 
model by using the growth function obtained in the previ- 
ous section and the initial condition for the curvature per- 
turbation and the scalar sp ectral index, = 2. 41 x 10 -9 , 
n = 0.96, given by WMAP l|Komatsu et alj|2009h . 

The 2D Fourier transform of the projected Compton 
2/-parameter is given by 

. sin(£x/£ s ) 



4-Kr B 



&xx y 3 d(x)- 



(30) 



where yzd is the radial profile of the Compton t/-parameter, 

(31) 



U3d(x) = — n e (x)k B T fi (x), 

with ctt being the Thomson cross section and m e the elec- 
tron mass. l s is the angular wavenumber corresponding to r s , 
£ s — Da/ts, with Da being the angular diameter distance. 

For the electron density pro file n e and the temperatur e 
profile T e , we use the result of iKomatsu fc Seliakl (2002), 
which is based on the NFW dark matter density profile 
l|Navarro et al.ll 19971 ). The NFW dark matter density profile 
is given by 



Pdm{x) — 



(32) 



where x := r/r B where r B is a scale radius, and p s is a scale 
density. The scale radius r s is related to the virial radius by 
the concentration parameter c as 



3 (M,z) 



r(M,; 



c(M,z) 



(33) 



We us e the concentration parameter of IKomatsu fc Seliakl 
(2003), 



3 (M, 



10 



1 



M 



M„(0) 



(34) 



where M*(0) is a solution to a(M,0) = S c at the redshift 

2 = 0. 

We adopt the spherical collapsed model to obtain the 
virial radius, 



3M 



AtyA v (z)p(z)' 



(35) 



where A v (z) is the virialised overdensity at z. Modified grav- 
ity chan ges the virialised ove rdensity, which we compute 
following ISchmidt et al.l (|2009l ). For example, the modified 
virial overdensity at z — is A„ = 270 for (a,r c Ho) = 
(0.1,0.3) and A„ = 300 for (a,r c H ) = (0.1,0.4), while 



A„ = 370 in the ACDM model. Thus, we find that the actual 
density in a virialised halo, A v (z)p(z), is lower in modified 
gravity than in the ACDM model. 

In order to o btain the profiles of the electron density 
and temperature. IKomatsu fc Seli ak (2002) assumed three 
things: (i) the electron gas is in hydrostatic equilibrium in 
the dark matter potential; (ii) the electron gas density fol- 
lows the dark matter density in the outer part of the halo; 
(iii) the equation of state of the electron gas is polytropic, 
P e oc p2, where P e , pe and 7 are the electron gas pres- 
sure, the gas density and the polytropic index, respectively. 
Under these assumptions, the electron number density and 
temperature profiles are simply given by 



n ec F(x), 



(36) 
(37) 



Here, n ec and T ec are the electron density and tempera- 
ture at the centre, respectively, and the dimensionless profile 
F(x) is written as 



F(x) = \ l - A 



1 - 



ln(l + c) 



1/(7-1) 



where the coefficient A is defined as 



A := 3??" 1 



7-1 



7 



ln(l + c) 



(38) 



(39) 



1 + c 

The central electron density n ec and the temperature T, 
are given by 

n ec = 3.01 



lO 14 M y \1 Mpc J 
1 -1 



(1 + c) 2 
T ec = 0.88 Vc 



'M/(W 14 h- 1 M @ ) 



F _1 (c) cm" 3 , (40) 
keV. (41) 



rvir/(l/i- x Mpc) 

For 7 and r] c , IKomatsu fc Seliakl (|2002l ) provided the follow- 
ing useful fitting formulas: 

7 = 1.137 + 8.94 x 10~ 2 ln(c/5) -3.68 x 10" 3 (c-5), (42) 
77c = 2.235 + 0.202(c - 5) - 1.16 x 10" 3 (c - 5) 2 . (43) 

Figure [6] shows the SZ power spectrum in the modi- 
fied gravity model. We take three sets of the model param- 
eters: (a,r c H ) = (0.1,0.2) , (0.1,0.3), and (0.1,0.4), giving 
as = 1.2, 1.0 and 0.9, respectively (The initial density power 
spectra are the same for all of these parameter sets). For 
comparison, we plot in Fig. [S]the SZ power spectra in the 
ACDM models with erg = 0.77 (our fiducial model) and with 
as = 1.0. One sees that the amplitude of the peak in modi- 
fied gravity is lower than that in the ACDM model with the 
same as- This is because the actual density in a virialised 
halo is lower in modified gravity than in the ACDM model. 

Moreover, modified gravity shows the damping of the 
SZ power spectrum on small scales and the amplification on 
large scales, because the contribution from halos at high red- 
shifts in modified gravity is smaller than that in the ACDM 
model. The redshift contribution for different I is shown in 
Fig. [7] The growth of matter density fluctuations in the mod- 
ified gravity model is rapid at late times compared with the 
growth in the ACDM model with the same as- As a result, 
the formation of halos delays, leading to the large contri- 
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Figure 6. SZ angular power spectra for different modified grav- 
ity parameters. The dotted, solid, and dashed lines represent SZ 
power spectra for (a, r c H ) = (0.1,0.2), (0.1,0.3), and (0.1,0.4), 
respectively. The SZ angular power spectra for as = 1.0 and 
erg = 0.77 in the ACDM model are shown as the thin solid and 
thin dashed lines, respectively. For references, we plot the primor- 
dial CMB temperature angular power spectrum in our fiducial 
ACDM model and the ACBAR data. 

bution from low redshifts and small contribution from high 
redshifts. 

The SZ power spectrum is constrained by the CMB ob- 
servations on small scales. In the ACDM model the ACBAR 
data at I ~ 2500 gives the c onstraint on the SZ power spec- 
trum in terms of erg: a 8 < 1 (|Kuo et al.ll2007h . The peak am- 
plitude of the SZ power spectrum is lower in modified grav- 
ity than in the ACDM model with the same erg. This result 
then yields the constraint from the ACBAR data as < 1.2 
in modified gravity. Figure [S] shows erg as a function of the 
model parameters a and r c Ho. In the limit r c — * oo, as re- 
duces to the value in our fiducial model, erg = 0.77. From 
the ACBAR experiment, the parameter region below the 
as = 1.2 line in Fig. [8] is excluded. 

Stronger constraints can be obtained from the small 
angular scale CMB me asurements by the QUaD experiment 
jFriedman et al.ll2009h . The QUaD telescope measured the 
CMB temperature anisotropy in the multipole range 2000 < 
£ < 3000 at 150 GHz. The QUaD team reported no strong 
evidence of the SZ effect and the result is consistent with 
the WMAP data which yields erg ~ 0.8. Therefore, using the 
result of the QUaD experiment to constrain as, the modified 
gravity parameters are more strongly restricted than in the 
case of the ACBAR data. The constraint from the QUaD 
data gives as ~ 0.9 in modified gravity. 



5 CONCLUSIONS 

In this paper, we have explored observational consequences 
of structure formation in modified gravity. The model 
considered is a phenomenological extension of the DGP 
braneworld, and modification to the standard ACDM model 
is characterized by the additional term ±H 2a /rt in the 
Friedmann equation. In the case of a — 1/2, the term arises 
from the five-dimensional effect, but in this work a and r c 
were assumed to be free parameters. 



I = 1000 

£ = 5000 / V 

' f w* *>:' '■• 

I = 9nnnn i i _ -\ . . 




0.1 1 

z 

Figure 7. Distribution of the redshift contribution of the SZ 
angular power spectrum for different i modes. We set the mod- 
ified gravity parameters a = 0.1 and r c Hg = 0.3. The solid, 
the dashed, and the dotted lines represent the distributions for 
£ = 1000, £ = 5000, and £ = 20000, respectively. For comparison, 
we plot the distributions for the ACDM model with erg = 1.0 as 
thin lines. 




Figure 8. The amplitude erg as a function of modified gravity 
parameter a and r c . The dashed-dotted-dotted, dashed, solid, 
dotted, dashed-dotted lines represent erg = 0.8, 0.9, 1.0, 1.1 and 
1.2. 

First, we have studied the spherical collapse model of 
nonlinear structure formation in modified gravity. It was 
found that change in the growth function is relatively large, 
but the linear density contrast for spherical collapse un- 
dergoes very small modification. For the self-accelerating 
branch, the growth of perturbations is suppressed compared 
to the ACDM model, which conf irms qualitatively the re- 
sult of ISchaefer fc Kovamal (|2008). For the normal branch, 
the growth of perturbations is enhanced compared to the 
ACDM model. 

Focusing on the normal branch, we then investigated 
the effect of modified gravity on the SZ angular power spec- 
trum. The enhanced growth function in the normal branch 
results in the amplification of erg. However, the modification 
to the SZ power spectrum is rather nontrivial. The peak 
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amplitude of the SZ angular power spectrum is lower in 
modified gravity than in the ACDM model with the same 
as, because halos are virialised at lower density in modified 
gravity than in the ACDM model. In addition to this, mod- 
ified gravity shows the damping of the SZ power spectrum 
on small scales. 

We confronted the modified SZ spectrum with CMB 
observations on small scales. Observational constraints can 
be read off from Fig. [5] The ACBAR experiment gives the 
constraint on the SZ power spectrum which translates to 
as = 1.2 in modified gravity. From this we put constraints to 
the parameters of the modified gravity model: for example, 
we have r c Ho > 0.4 for a = 0.5 and r c Ho > 0.2 for a — 0.1. 
The QUaD data gives severer constraints on the parameters: 
as = 0.9, leading to r c H > 0.4 for a = 0.1. 




0.2 0.4 0.6 0.8 1 

a 

Figure Al. Comparison with the growth functions obtained by 
the different approaches. 
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APPENDIX A: COMPARISON WITH 
DIFFERENT APPROACHES 

In the phenomenologically extended version of the DGP 
model, different assumptions can be made in computing the 
growth of perturbations. Instead of using Birkhoff 's theo- 
rem, we may assume that the behaviour of pert urbations 
is gov erned by a scalar-tensor theory. Following iKovamal 
(2006), we now compare the linear evolution of perturba- 
tions obtained by these two approaches. If one does not as- 
sume Birkhoff 's theorem but employs a scalar-tensor theory 
as an effective theor y for perturbat ions, the evolution of 8 
will be described by l|Kovamal 120061 ') 



.i + 2HS = 4ttG ( 1 + — ) pS, 



where 



1 + 



(rcH) 



2(1-0,) 



1 + 



2 (l-a)H 



H 2 



(Al) 



(A2) 



Note that this is the expression for the normal branch. 
We denote by D+(a) the growth function calculated from 
this equation. The difference between D\ and D+ obtained 
in the main text, (D+ — D+)/(D* + + D+), is presented in 
Fig. IA1I As can be seen, there is only a few percent dif- 
ference between the two approaches concerning the linear 
growth of perturbations. Developing nonlinear theory with 
the violation of Birkhoff 's theorem taken into account is be- 
yond the scope of the present paper. 
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